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ABSTRACT 

We present a multi wavelength study of the interacting starburst galaxy NGC 7771, 
including new optical and ultra-violet spectra and a previously unpublished soft X-ray 
ROSAT image and spectrum. The far-infrared, radio, and X-ray fluxes suggest that a 
massive burst of star-formation is currently in progress but the small equivalent width 
of the Balmer emission lines (equivalent width Ha « 100 A), the weak UV flux, the 
low abundance of ionised oxygen, and the shape of the optical spectrum lead us to 
conclude that there are few O stars. This might normally suggest that star-formation 
has ceased but the galaxy’s barred gravitational potential and large gas reserves imply 
that this should not be so, and we therefore consider other explanations. We argue that 
the observations cannot be due to effects of geometry, density bounded nebulae, or dust 
within the nebulae, and conclude that a truncated IMF is required. The dwarf galaxy 
NGC 7770 appears to be in the initial stages of a merger with NGC 7771, and the 
resulting tidal perturbations may have induced the apparent two-armed spiral pattern, 
and driven a substantial fraction of the disk gas inwards. The presence of a bulge in 
NGC 7771 may be moderating the starburst so that, while still occuring on a large scale 
with a supernova rate of 0.8-1 yr _1 , it is less violent and the IMF has a relatively low 
upper mass limit. We find that there is a cluster of stars obscuring part of the starburst 
region, and we offer an explanation of its origin. 

Key words: Galaxies: starburst - Galaxies: individual: NGC 7771 


1 INTRODUCTION 

The term ‘starburst’ was originally coined by Weedman et 
al. (1981) to indicate qualitatively that a galaxy is undergo- 
ing a period of intense star formation. It still lacks a precise 
definition and provides a blanket covering for all regions that 
have an above normal star formation rate. In particular for 
interacting galaxies, it includes star formation enhanced by 
any process ranging from a mild perturbation to a collision. 
When the data available on a certain system is limited and 
the initial mass function (IMF) must be guessed at, often the 
assumption is that the starburst has a Salpeter IMF with 
masses in the range 0. 1-100 M©. While this may be a reason- 
able starting point for analysis of starbursts, it has not been 
established that it can be applied directly to all star-forming 
regions. Indeed there is evidence to the contrary. 

There has been some effort to constrain the IMF mass 
limit and Chariot et al. (1993) have investigated the evolu- 
tion of starbursts with high lower mass cutofTs. Scalo (1989) 
reviewed the evidence for top-heavy IMFs resulting in a de- 


ficiency of low-mass stars in starbursts. He concluded that 
most starburst galaxies do not require a top-heavy IMF, 
if the upper mass limit (M u ) is large enough (M u > 80- 
100M©). In a significant fraction of galaxies this was not 
the case and so to avoid excessive masses of newly formed 
stars the lower mass limit (Mi) must be higher, perhaps as 
much as Mj = 3-1 OM© as was found by Rieke et al. (1980) 
for the starburst in M82. Which of these situations applies in 
any particular case has bearing on the masses involved and 
the efficiency with which stars are formed, as well as the 
question of whether perturbations and heating associated 
with enhanced star formation rates inhibits the formation 
of low-mass stars. More recently Doyon et al. (1992) sug- 
gested that the ratio Hel 2.058^m/Br7 might be a reliable 
indicator of M u for ratios less than ~0.5, corresponding to 
M u < 40M©, although Shields (1993) pointed out that the 
neutral helium to hydrogen fraction decreases with increas- 
ing effective temperature. 

Evidence for truncated IMFs in starburst galaxies is 
very limited at present, due partly to the number of factors 
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which must be considered, and also because many starbursts 
may indeed have IMFs that are not significantly truncated. 
Perhaps the best documented case of a truncated IMF is in 
M82, for which the upper mass cut-off has been found to 
be about 30 Mq (Rieke et al. 1980, Puxley et al. 1989, and 
Doyon et al. 1992); more recent work on this galaxy by Rieke 
et al. (1993) using updated models has shown that while low 
mass star formation appears to be suppressed, the abun- 
dance of higher mass stars is consistent with that found in 
the solar neighbourhood. The debate on this galaxy contin- 
ues and, by measuring mid-infrared line ratios, Achtermann 
& Lacy (1995) place an upper limit on the effective temper- 
ature of the hottest stars as ~35,000 K, equivalent to stars 
of ~20 M©. The argument for unusual IMFs been further 
supported by Gehrz et al. (1983), Wright et al. (1988), and 
Doyon et al. (1992) who have found other galaxies where the 
mass function is severely truncated at both ends. 

In this paper we consider the case of NGC 7771, an 
SB(s)a galaxy interacting with NGC 7770, at a redshift of 
4345 km s _1 (Soifer et al. 1987) corresponding to a distance 
of 87 Mpc ( Ho — 50 km s" 1 Mpc -1 ). This galaxy has an 
infrared luminosity of log L/Lq = 11.6 ranking it amongst 
some of the most luminous infrared galaxies known. The 
stellar origin for this power output and the galaxy’s rela- 
tive closeness with respect to other comparably luminous 
objects, makes it an ideal candidate for the study of star- 
formation on a massive scale. Our observations reveal that 
the central region of NGC 7771 has unusual characteristics 
which could be explained in a variety of ways, one of which 
requires a truncated IMF. 

This paper is organised as follows: In Section 2 we 
present and discuss the optical spectrum, in Section 3 the 
UV spectrum, and in Section 4 the X-ray emission. These 
are put into context with each other and discussed with ref- 
erence to detailed starburst models in Section 5. Finally, in 
Section 6 we present our conclusions. 

2 OPTICAL SPECTRUM 
2.1 Observations 

The optical spectrum was obtained using the Intermediate 
Dispersion Spectrograph (IDS) on the 2.5-m Isaac Newton 
Telescope (INT) on November 12, 1993 with a 900s inte- 
gration. The configuration was a 1.5 arcsec wide long-slit 
centered on the optical peak and oriented at position angle 
59° , shown schematically in Figure 1. The spatial resolu- 
tion along the slit was 0.6 arcsec/pixel although the PSF of 
the standard star showed that the seeing was about 1.5 arc- 
sec. The spectral resolution of R ~ 1000 was matched to 
the 3.05 A/pixel, yielding an effective resolution of about 
6 A (FWHM) and covering the range approximately 3600- 
6900 A. 

Data reduction was performed with I RAF, following 
standard procedures. After bias subtraction and flatfielding, 
the pixels affected by cosmic rays were rejected and inter- 
polated over. The background was fitted using regions away 
from the bright galaxy, and subtracted. Wavelength cali- 
bration was carried out using a Cu-Ar arc lamp and then 
spectra were extracted, corrected for atmospheric extinc- 
tion, and flux calibrated by a standard star observed at a 
similar air mass. 



Figure 1. Contour plot of NGC 7771 and NGC 7770 in the opti- 
cal from the Digitised Sky Survey showing the orientation of the 
1.5 arcsec wide long-slit. The shaded region of the slit shows the 
section that has been mapped in Figure 2; the two bars parallel 
to the slit indicate the two regions discussed in Section 2.3. 


2.2 Spatial Mapping 

In order to make the most use of the spatial information 
that a long-slit affords, the central 19.5 arcsec was divided 
into a series of 1.5 arcsec (630 pc) apertures. For each of 
the resulting spectra measurements were made of the total 
integrated flux of the spectrum, the slope of the continuum, 
the Ha line flux and equivalent width (EWh q ), and where 
possible the H/? line flux. 

Figure 2 shows in (a) the total brightness, computed as 
simply the total flux detected in each aperture normalised 
to 10 for the central one, and in (d) the continuum colour, 
i.e. the ratio of the continuum flux density at Ha to that at 
H/? with the dotted line representing a flat continuum. Both 
quantities show clearly that along the slit, the region within 
±4 arcsec of the peak, about 3 kpc across, is brighter and 
bluer, indicative of a starburst nucleus. The red outer regions 
are consistent with an old evolved stellar population similar 
to that found in galactic bulges. The Ha flux and equiv- 
alent width, in (b) and (e) respectively, show immediately 
that this is not the complete picture. In both diagrams the 
optical brightness has been over-plotted with a dashed line 
and arbitrary scale for comparison. To explain why these do 
not peak at the same place, we first consider the northeast 
side (negative offset). The Ha flux increases dramatically at 
the same point as the optical brightness, which is also where 
the continuum becomes much bluer. This suggests that the 
nucleus hosts a cluster of hot young OB stars which are pro- 
ducing a large ionising flux. In such a situation we might 
expect EWho to change more slowly, which is the case, 
as further from the nucleus the continuum becomes more 
dominated by older redder stars. However, to the southwest 
(positive offset) the Ha flux decreases abruptly. That the 
EWhq also suffers a marked cut-off suggests that the stellar 
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Figure 2. Mapped region of the long-slit showing the variation along the slit. Negative offsets are to the northeast and positive offsets 
southwest. The individual frames (described in the text in more detail) are: (a) total optical brightness, (b) observed H a flux, (c) ratio 
Ha/H/3, (d) continuum colour, and (e) Ha equivalent width. 


population changes suddenly from hot OB stars to cooler 
non-ionising ones, but this cannot be correct because the 
continuum colour shows that the stars are still very blue. 
To understand it further we need to know how the extinc- 
tion varies across the nucleus. Since this can be calculated 
from the Balmer decrement, we have shown the ratio Ha/H 0 
in (c), although it could only be measured for the 5 aper- 
tures in which H/? was strong enough to detect. The Ha 
flux has been over-plotted (dotted line and arbitrary scale) 
to show that it is strongest where the extinction is lowest 
(Av ~ 2 mag), and that the cut-off occurs where the extinc- 
tion increases rapidly by 2-3 magnitudes. The data plotted 
are maximum ratios as there is considerable H/? absorption, 
especially to the southwest, which could not always be ac- 
curately deblended from the emission line flux. It is notable 
that the dereddened Ha flux is approximately constant over 
the region where it can be calculated, suggesting that the 
starburst does indeed cover the entire region indicated by 
the bluest continuum colours. There might then be either a 
dust lane or a cluster of non-ionising stars obscuring part of 
this, and causing the extra extinction. The effect of a dust 
lane would be to redden the continuum and reduce the de- 
tected flux. Since neither of these happen, and since also 
there is no evidence for a dust lane on the optical image, we 
conclude that the obscuration is due to a cluster of stars. In 
order to investigate these two regions with a higher signal- to- 
noise ratio, we extracted two further non-overlapping spec- 
tra, both 2.4 arcsec wide. The first, denoted region 1, was 
offset to the northeast by 3 arcsec, to look at the starburst 


where it was least obscured; the second, denoted region 2, 
was centered on the optical peak, corresponding to the ob- 
scuring cluster and high extinction to the starburst. These 
two regions are indicated in Figure 1 by the two thick lines 
drawn parallel to, and within, the slit. 

Although the obscuration clearly varies considerably 
across the mapped section of the nucleus, the average ex- 
tinction is a useful quantity which we estimate to be Av = 
2.3 mag, noting that it will be biassed towards lower val- 
ues that are more easily measured. We have shown that the 
starburst, although partly obscured, extends with approxi- 
mately constant surface brightness across the whole 7.5 arc- 
sec nucleus at the orientation of the slit. In this limited re- 
gion, we calculate an extinction corrected surface brightness 
for the Ha emission of 4.0 x 10“ 16 W m -2 kpc -2 . 

2.3 Two Nuclear Regions 

The spectra for both regions, smoothed with a Gaussian fil- 
ter of FWHM ~6 A are presented in Figure 3 on the same 
scale to allow direct comparison. The most obvious differ- 
ences are that Region 1 is less reddened and Region 2 has 
greater Balmer absorption as well as a bluer continuum. 

From the Balmer emission line decrement we find that 
in Region 1 Av = 1.9 mag, while in Region 2 Av = 3.0 mag. 
The extra obscuration is caused by the star cluster that is 
also responsible for the much enhanced Balmer absorption. 
The spectrum of this cluster is discussed in more detail later 
(see Section 2.4). The strengths of the detectable emission 
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Figure 3. Spectra of the starburst where it is least obscured, Region 1, and where the obscuring cluster is covering it, Region 2. They 
are plotted on the same scale to emphasise the differences in both lines and continuum. 


Table 1. Emission Lines 


emission 

line 

F 1 

r oba 

Region 1 
F h2 
r corr 

EW (A) 

F 1 

r obi 

Region 2 

p 1 > 3 
1 corr 

EW (A) 

[Sll] 

A6731 

5.1 

20 

8.7 

2.3 

20 

2.6 

[SII] 

A6716 

5.6 

23 

9.8 

3.3 

29 

3.7 

[Nil] 

A6583 

23.8 

99 

42.0 

11.9 

108 

13.2 

Ha 

A6563 

60.3 

253 

106.5 

28.6 

262 

4 31 .8 

[Nn] 

A6548 

8.1 

34 

14.4 

5.1 

47 

5.6 

[Ol] 

A6300 

0.8 

4 

1.4 

0.8 

8 

0.9 

[Hel] 

A5876 

1.0 

5 

1.6 

0.2 

3 

0.2 

[Olll] 

A5007 

1.0 

7 

1.5 

0.7 

15 

0.7 

Hp 

A4861 

10.8 

87 

16.2 

3.6 

89 

3.2 

H-7 

A4340 

2.7 

30 

4.0 

— 

— 

— 

[On] 

A3727 

3.0 

52 

5.9 

— 

— 

— 


1 observed and reddening corrected fluxes xlO“ 18 Wm~ 2 in a 2.4 arcsec x 1.5 arcsec aperture. 
la errors, not including deblending, are typically ±0.3 in the same units. 

2 1*7 errors typically ±2 at H (3 and ±1 at Ha. 

3 la errors typically ±6 at H j3 and ±2.5 at Ha. 

4 Veilleux et al. (1995) find that EWn a = 26 A in a 2 arcsec X 2 arcsec aperture centred on the optical peak. 


lines before and after appropriate dereddening, as well as 
their equivalent widths, are presented in Table 1. 

Kim et al. (1995) find Ha/H /? = 21.4 in a 2 arcsec 
x 2 arcsec aperture centred on the nucleus, considerably 
larger than our value of 8.1 for Region 2, due in part to our 
measuring a larger equivalent width, 4.6 A, for the H/? ab 
sorption. We do, however, measure a similar Nal absorption 
of 3.3 A in the same region. 

A simple test can be performed to determine whether 
the emitting gas is photoionised by OB stars or a power- 


law continuum source. We employ the diagnostic diagrams 
described in Fruscione h Griffiths (1991), which compare 
the dereddened ratios shown in Table 2. Plotting the line 
ratios on the diagrams shows clearly that both regions are 
starbursts. Since the Ha surface brightness is the same and 
they are in close proximity with no obvious break between 
them, it seems reasonable to conclude that they are in- 
fact part of the same starburst. The difference in reddening 
(AAv =1.1 mag) is entirely attributable to the obscuring 
cluster. Storchi-Bergmann et al. (1995) found for a sample of 
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Table 2. Emission line ratios 


ratio 

Region 1 

Region 2 

log [Om]A5007/H/? 

-1.08 

-0.77 

log [Nll]A6583/Ho 

-0.41 

-0.38 

log [Sll](A6716 + A6731)/H« 

-0.77 

-0.73 

log [Ol]A6300/Ho 

-1.85 

-1.50 


Note: line ratios are corrected for reddening. 


starbursts that log [Om](A4959 + A5007)/H/? = -O.lOlJgJ 
(not corrected for extinction). Our values (which are the 
same within their la errors) lie at the bottom end of this 
range, and are consistent with an effective temperature for 
the radiation field of the starburst of 35-40 xlO 3 K, assum- 
ing solar metallicity (Evans &; Dopita 1985). 

An upper limit to the velocity dispersion implied by 
the widths of the lines can be estimated by comparing their 
FWHM to the instrumental FWHM of 287±8 km s -1 (mea- 
sured from the arc lines). We find that for Region 1 the 
FWHM of the six strongest lines is 291 ± 9 km s -1 , showing 
that they are unresolved; and for Region 2 the FWHM for 
the same lines is 370 ± 12 km s _1 . The limit on the velocity 
dispersion is therefore <200 km s _1 , in line with what we 
would expect for a nuclear starburst; we suggest a reason 
for the differing line widths in Section 5.4. The lack of any 
broad component suggests there is no outflow such as from 
supernovae induced winds (see Section 5.2). The heliocentric 
velocity of the lines is 4790 ±80 km s~ J . 

The ratio of the sulphur lines in the doublet can be used 
as a density diagnostic. We find that [Sll] =1. 1-1.4, 

typical of densities 100-500 cm -3 , consistent with those val- 
ues found in other extragalactic Hll regions (French 1980). 

2.4 The Obscuring Cluster 

Since the emission lines in the two regions are produced 
by the same starburst, we can redden the spectrum from 
Region 1 using AAv = 1.1 mag so that the emission lines 
correspond to those in Region 2 and then subtract it without 
further scaling. The residual will be the spectrum of the stars 
and dust that is obscuring the starburst in Region 2, and is 
shown in Figure 4 smoothed with a Gaussian of FWHM 9 A. 

The deep Balmer absorption is clearly shown. There 
also appears to be some residual Ha due to the difference 
in line-widths between the two regions; the overall excess 
emission in the line is very small. The continuum is very 
blue: the ratio of the continuum at Ha to that at Hy? is 
0.68, more extreme than any in Figure 2(d). There can be 
very little foreground extinction since this would make the 
intrinsic spectrum even bluer, which would lead us to expect 
some of the stars contributing to it would be hot enough to 
create an Hll region. In such a situation the Ha flux and 
equivalent width would not drop as abruptly as is seen. 

We have used the stellar library of Silva & Cornell 
(1992) to synthesise this spectrum. By considering partic- 
ular features of the spectrum it is possible to constrain 
the range of populations in the model. For example, the 
Balmer decrement (3700-4000 A) is typical of A and B stars, 
whereas the absorption bands of Mgl (5150-5200 A) and 
Nal (5880-5900 A) are indicative of K and M stars; and the 
Cal (4230 A) and Call (3900-4000 A) absorption is most ap- 


parent in F stars. Lastly, the overall slope of the continuum 
must be matched and puts further limitations on the possi- 
ble stellar classes present. Although this does not provide an 
entirely unique model, it does show broadly the fractional 
contributions from the various stellar types. The spectra in 
the stellar library are normalised at 5450 A and it is at this 
wavelength that we find significant contributions from B, A, 
F, and G stars, but also enough optically faint K and M 
stars that these will dominate the mass. A small amount of 
reddening (Av = 0.3 mag) is also required. The H/?, H7, 
and Htf absorption features indicate that the redshift of the 
object is 4800 ± 240 km s"" 1 , the same as the redshift of the 
emission lines in Section 2.3. 


3 UV SPECTRUM 

A UV spectrum was obtained with the International Ul- 
traviolet Explorer (IUE) satellite on November 29 and 30 
1994 with the long-wavelength (LWP) and short- wavelength 
(SWP) cameras respectively, in the large (10 arcsec x 
20 arcsec) aperture. The resulting spectrum, extending from 
3200 A to 1200 A, was extremely faint and so we were unable 
to distinguish or measure any line fluxes. In order to increase 
the signal-to-noise ratio we binned the data in 200 A sec- 
tions, 3 each for the LWP and SWP cameras. Because of this 
we have not shown the spectrum here explicitly although it 
is included in Figure 8 (see Section 5.1). 

Kinney et al. (1993) looked at a sample of starbursts 
for which p = — 1.25± 0.45 (F\ = f 0 X 0 ). NGC 7771 has p = 
0.56 ±0.18 (weighted towards the fluxes with higher SNR), 
considerably redder than even the reddest in the sample of 
Kinney et al. (NGC 7552 with P = 0.26 ± 0.14). One of the 
criterion for objects to be included in the sample was based 
on signal-to-noise ratio and hence the sample is mostly of 
galaxies with low reddening and strong UV flux. The UV 
spectra of objects in the sample are consistent with those 
of O and B stars. Naively then, NGC 7771 would appear 
to be a heavily reddened source. However, Kinney et al. 
discuss the effects of reddening on the UV spectrum. They 
claim that the optical depth is unlikely to be near unity. 
Rather the dust will be patchy, and optically thick enough 
that either the UV radiation has a completely free path or 
encounters dust and is totally absorbed. Thus the effect of 
dust is not to alter the shape of the spectrum, but simply to 
reduce the observed flux. This is supported by their inability 
to reproduce a typical blue spectrum by dereddening one of 
the red spectra, even with a variety of extinction laws. We 
conclude that the shape of the spectrum in NGC 7771 is not 
due to reddened flux from O and B stars, but is from A or 
B stars which it more closely resembles. 


4 X-RAY DATA 

Fabbiano (1989) suggests that the typical soft X-ray lu- 
minosity of starburst nuclei is on the order of 10 32 W. 
NGC 7771 is unusual in having an X-ray luminosity of more 
than 10 34 W (see below), and so we have included a section 
on the observations and possible origins of the emission. 

NGC 7771 has been observed in soft X-rays with both 
the Einstein Observatory and ROSAT, but no hard X-ray 
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Figure 4. Spectra of the obscuring cluster, obtained by subtracting the suitably reddened spectrum of Region 1 from that of Region 2. 


emission was detected by the HEAO-1 mission; for example 
the data presented in Rephaeli et al. (1995) have less than 
1<7 significance. All X-ray data includes contributions from 
NGC 7770 as well as NGC 7771 and we have not attempted 
to distinguish between them. 

We have recalibrated the soft X-ray flux detected by the 
Einstein IPC device (0.16-3.5 keV) from the CDROM im- 
ages. In a region with radius 210 arcsec we measure 100.2±10 
counts, agreeing well with the 104.7 measured by Fabbiano 
et al. (1992) for the same aperture. The hardness ratio, that 
is the relative contributions of ‘hard’ (0.81-3.5 keV) and 
‘soft* (0.16-0.81 keV) counts, for the galaxy is 0.36 ± 0.13. 
Together with the Galactic absorbing column of hydrogen, 
Nh = 3.96 x 10 20 cm" 2 , it can be used to estimate the flux 
using an appropriate model. We choose the Raymond-Smith 
thermal plasma model with cosmic abundances and a tem- 
perature kT = 1.7 keV (from the hardness ratio) to give a 
luminosity of 5.4 ± 0.5 x 10 34 W. 

The galaxy was also observed by ROSAT on two occa- 
sions with the PSPC device: on December 16 1992 with a 
live-time of 7655 s, the count rate was 0.0151 ± 0.0020 s -1 
in an aperture of radius 210 arcsec; and on July 10- 
13 1992 with a live-time of 4150 s, the count rate was 
0.0181 ±0.0043 s" 1 for the same aperture. 

The ROSAT PSPC X-ray images from both observa- 
tions were aligned (with reference to several prominent ob- 
jects in the field, and requiring sub- pixel shifts) and aver- 
aged, weighting each frame with its exposure time. The re- 
sulting image is shown in Figure 5 as contours (linear scale) 
superimposed on a grey-scale optical image. The 30-arcsec 
resolution of this image is better than that of the Einstein 
IPC (75 arcsec) and appears to show that, although the 
strongest emission is from the galactic centre, it extends 
throughout the galaxy’s disk and possibly beyond. The emis- 
sion to the south of the galaxy and the ridge which appears 
to join NGC 7770 may be effects caused by the PSF of the 
PSPC which is quite large with respect to the angular size 
of the galaxy (it is represented by the size of the scale bar). 
Emission from near the edge of the disk may appear to origi- 
nate from outside it. Higher resolution imaging, such as with 
the HRI instrument, is needed to resolve this point. 

The PSPC instrument has 256 channels which provide 
a certain amount of spectral coverage (0.1-2. 4 keV). For 
the latter observation above we extracted this spectrum us- 



Figure 5. Contours showing the X-ray emission detected by the 
ROSAT PSPC, centered on a grcy-scale optical image taken from 
the Digitised Sky Survey. The PSPC has a FWHM resolution of 
30 arcsec, as shown by the scale bar. 


ing ASTERIX and corrected both source and background 
counts for vignetting and dead time. The background was 
subtracted using XSPEC, and channels 20-200 of the result- 
ing spectrum were modelled. Due to the small total counts 
(only about 100 photons were detected) the spectral analysis 
is restricted and we do not quote formal errors. Instead the 
‘best-fit’ models should be taken as in indication of the possi- 
ble fits rather than as unique solutions. Two models, statisti- 
cally indestinguishable by minimum \ 2 (reduced x 2 = 0*4), 
are shown in Figures 6 and 7 together with the data which 
has been binned to a level of 2<r or a maximum of 16 chan- 
nels per bin. Figure 6 shows a power-law fit for which the 
best fitting parameters are an absorption column density 
JVh = 1.8 x 10 21 cm" 2 , and a photon index T = 3.5. Fig- 
ure 7 shows in (a) a fit with an absorption column density 
Ah = 4.5 x 10 2 ° cm -2 (similar to the galactic absorbing 
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Figure 6. ROSA T PSPC X-ray spectrum (binned to a level of 2 a 
or a maximum of 16 channels) shown with the best-fit power- law 
model with absorption. 



0.4 0.6 0.6 t 2 

channel energy (keV) 

Figure 7. ROSAT PSPC X-ray spectrum (binned to a level of 
2 a or a maximum of 16 channels) shown with best-fit Raymond- 
Smith hot plasma model with absorption, (a) shows the model, 
and (b) shows the residuals which indicate line emission at 0.65 
keV. 

column of Nu = 3.96 x 10 20 cm -2 ) and a Raymond-Smith 
model with cosmic abundances and a plasma temperature 
given by kT = 0.6 keV (in good agreement with the canoni- 
cal model described in Serlemitsos et al. (1996) for a sample 
of LINERs and starburst galaxies). In (b) we show the resid- 
uals which are consistent with line emission from Ovm at 
0.65 keV and the Fe-L complex at 0. 8-1.0 keV (Serlemit- 
sos et al. 1996). Using these models to calibrate the ROSAT 
counts we find a luminosity for the power-law of 1.6 xlO 34 W, 
and for the plasma model 1.2 x 10 34 W. The difference be- 
tween these and the luminosity calculated from the Einstein 
data may be due partly to the uncertainty of the model 
parameters. 

We consider the merits of these two models and their 
relevance for starbursts. A power-law model is used for 
analysing Seyfert galaxies, when the photon index is typ- 
ically r = 2.3 ± 0.7 (Boiler et al. 1992). The emission line 
ratios and extended nature of the X-ray flux show that this 
is clearly not the case here. The X-ray emission in star- 


burst galaxies may be powered by massive X-ray binaries, 
for which we would expect 0.8 < T < 1.5 (Nagase 1989). Na- 
gase also points out that the associated absorbing column 
densities are typically very high (on the order of 10 23 cm -3 ). 
Our model predicts an absorbing column higher than galac- 
tic, but it is harder to reconcile the large photon index 
with that expected for X-ray binaries. We can estimate the 
star-formation rate required for this model by following the 
method of Ward (1988) and assuming the average luminos- 
ity of an X-ray binary to be 10 31 W. If we take the lifetime 
of such a binary to be 4 x 10 4 yr (Fabbiano et al. 1982) we 
find that a new binary is needed every 25 years. 

Alternatively, if the X-rays arise from supernovae, the 
Raymond-Smith model for hot plasma would be more appro- 
priate. We approach the problem by adopting an average X- 
ray luminosity of 10 30 W per supernova remnant, and a typ- 
ical lifetime of 5-10xl0 3 yr (Fabbiano et al. 1982, Rephaeli 
et al. 1995). The supernova rate is then 1.2-2. 4 yr -1 . 

Any realistic solution is likely to be a combination of 
both of these processes. In order to estimate the relative 
contributions, we can find the supernova rate that would be 
associated with the creation of a massive X-ray binary every 
25 years. If one X-ray binary occurs for every 500 OB stars 
(Fabbiano et al. 1982), then even if only 1/4 of these stars 
are supernova progenitors (assuming a Salpeter IMF in the 
range 3-100 M© for OB stars) we find that the X-ray emis- 
sion from supernovae must dominate over that from X-ray 
binaries. Although there will inevitably be a considerable 
scatter in the quantities we have used in our estimates as 
they depend strongly on the environment, we conclude that 
the parameters of the Raymond-Smith model, which is ap- 
plicable to supernova remnants, provide a more consistent 
solution to the interpretation of the X-ray emission. 


5 DISCUSSION 

5.1 Spectral Energy Distribution 

The spectral energy distribution is shown in Figure 8. Our 
calibration of the Einstein X-ray counts are represented by a 
bow-tie showing the uncertainty in the X-ray hardness. The 
optical and UV spectra described in this work are included 
as solid lines. The photometry in literature is as follows. 
The radio fluxes are from the Texas survey (365 MHz), Gre- 
gory Sz Condon (1991) (4.85 GHz), and Condon et al. (1990) 
(1.49 GHz). The spectral index between the 4.85 GHz and 
365 MHz fluxes (F u oc t'” 0 6 ), which both include emission 
from NGC 7770, is drawn as a dotted line. The 1.49 GHz 
flux of 124 mJy is integrated over NGC 7771 alone; includ- 
ing NGC 7770 would increase it by 16.5 mJy, not a sig- 
nificant amount as the line joining the first two fluxes also 
demonstrates. The sub-millimetre data is from Carico et al. 
(1992) (30 arcsec beam) and Hughes et al. (1997) (typically 
16.5 arcsec beam), and far-infrared from the Catalogue of 
Galaxies and Quasars in the IRAS Survey (Lonsdale et al. 
1988). The near-infrared data is from Carico et al. (1991) 
(5 arcsec aperture) and Eales et al. (1990) (1.4 arcsec and 
4.7 arcsec) The UBV photometry is from the 3 rd Reference 
Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991). 

It is immediately clear from Figure 8 that the far- 
infrared dominates the energy output with an enormous im- 
plied formation rate of massive stars. Some authors (e.g. 
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Figure 8. Spectral energy distribution of the galaxy. Circles denote fluxes measured through a large aperture, typically the whole galaxy; 
triangles denote smaller aperture (<20 arcsec) fluxes; the two solid lines show the new optical and UV spectra; the bow-tie represents 
the Einstein X-ray flux and the uncertainty of its hardness; the dotted line indicates a radio power-law. 


Sauvage fc Thuan 1992) have argued that a significant frac- 
tion of this luminosity is cirrus, and not related to such star 
formation at all. They claim that the decrease in Zfir/^Hcj 
with Hubble type for spiral galaxies depends on the rela- 
tive contribution from the cirrus component. Hughes et al. 
(1997) have shown that in NGC 7771 the coolest compo- 
nent of the dust re-radiation accounts for 3/4 of the total 
(8-1000 fim) infrared emission. If this is dominated by cir- 
rus, the nuclear starburst itself may be responsible for only 
1/4 of the infrared luminosity. However, due to the uncer- 
tainty and debate about this issue we have not tried to deter- 
mine the cirrus contribution. As the cirrus must be heated 
by starlight we can, by using the whole infrared luminosity 
from 8-1000/xm, measure the total star-formation rate in- 
cluding both the starburst and the continuous background 
rate for the galaxy. This will provide an estimate which may 
be more reliably compared to those derived from the total 
radio and total X-ray fluxes. 

A further complication arises because the resolution 
of the X-ray and far-infrared observations is insufficient to 
clearly distinguish NGC 7770 from NGC 7771. For the dis- 
cussion in the rest of this Section we have included both of 
these unless explicitly stated. 

An important quantity for starbursts is the supernova 
rate ( !/$*), reflecting the scale of the massive star formation, 
which we are able to estimate in several ways. We have al- 
ready seen that the X-ray flux implies i/ SN = 1-2, and we now 
look at that derived from the radio continuum. The radio 
spectral index between the 4.85 GHz and 365 MHz fluxes is 
a = 0.6 (F u oc consistent with the radio emission be- 

ing entirely non- thermal synchrotron radiation from shocks 
in supernovae. There will be a small component of thermal 
(free-free) emission in the ionised nebulae, but without a 


measurement of the total Ha flux we can only constrain t / s n 
with an upper limit. Condon and Yin (1990) discuss the 
supernova rate implied by the non-thermal radio flux and 
argue that 

, -U_. £nt(W Hz -1 ) 

Myr ) ~ L3 x 10 23 „-<*(GH z ) 

which yields i/ SN = 1.2 yr" 1 , or for NGC 7771 alone 1.0 yr” 1 . 

An alternative is to look at the total infrared flux (8- 
1000/zm), which can be found to within ~5% (ignoring un- 
certainties in the measurements) from the IRAS data us- 
ing the prescription of Perault (1989). For NGC 7771 and 
NGC 7770 combined it is 4.3 x 10 n .Z/©. To calculate the 
the star-formation rate (SFR) and supernova rate from this 
we will use the evolutionary synthesis models of Leitherer 
Sc Heckman (1995). Complications arise due to the metal- 
licity which can have significant affects, but in the absence 
of any evidence to the contrary we shall assume solar abun- 
dances throughout. The star- formation history is one of the 
most important parameters, particularly whether it is con- 
tinuous and if so how long it has been progressing. Although 
the starburst covers a region 3 kpc across, recent HST im- 
ages (eg Meurer et al. 1995) have shown that in many cases 
the star- formation occurs in smaller clusters. The overall ef- 
fect is for continuous star- formation over a larger region. 
The high supernova rates derived from the X-ray and radio 
fluxes imply that the star-formation has been progressing 
long enough for even 8-10 M© stars (which account for most 
of the supernovae) to evolve away from the Main Sequence. 
We therefore assume star-formation has been progressing for 
around 10® yrs. The models show that by this time, a sort of 
steady-state has been reached, and the parameters change 
very little with time. We find that for a Salpeter IMF in the 
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range 1-100 M®, SFR = 25 Mq yr*" 1 and i/ S n = 0.5 yr -1 
rather less than those derived from the radio continuum and 
X-ray flux. 

We can use i/sn to estimate the total Ha flux in 
NGC 7771 as 2-3 x 10~ 14 W m~ 2 . From our measurement of 
the extinction corrected Ha surface brightness in Section 2.2 
we find this would mean the starburst extends over an area 
on the order of 60 kpc 2 or around 20 x 20 arcsec, signifi- 
cantly more than is implied by the 7.5 arcsec extent of the 
line emission along the slit. Carico et al. (1988) have shown 
that the ratio for 10/xm flux within a 4.6 arcsec aperture to 
the IRAS 12 /xm flux (colour corrected to 10 fxm) is only 
0.17. Since the 10/xm emission is believed to be re-radiation 
from hot dust, it will trace the distribution of the OB stars 
heating it, and hence implies that the stars responsible are 
spread out. The radio morphology is also clearly extended 
along the major axis, although at a much weaker level than 
the nuclear emission. However, neither of these of these ob- 
servations can be reconciled with the implied extent of Ha 
emission, and so we consider the starburst in more detail. 

5.2 Equivalent Width of Ha 

The equivalent width (EW) of Ha emission is a useful di- 
agnostic of the starforming region. In particular it is inde- 
pendent of the extinction, and is less affected by stellar ab- 
sorption lines than is the H/? emission. Our observed value 
of EWhq ~ 100 A is somewhat lower than expected: a 
standard IMF with continuous star-formation will produce 
EWho > lOOA. We therefore discuss some possible causes 
for this discrepancy, beginning with the geometry. Recent 
HST images of starburst galaxies have shown that star- 
formation often occurs in many individual compact clusters. 
The effect of observing at low spatial resolution is for the 
presence of an underlying stellar continuum to reduce the 
equivalent widths of emission lines. However, in this case 
it cannot be significant as the underlying continuum con- 
tributes only 10-20 per cent of the total continuum in the 
region where we detect emission line flux (see Figure 2). The 
near-infrared spectrum of Smith et al. (1996) shows a low 
equivalent width for the Br 7 line, but this is expected since 
the contribution from an old population of stars is much 
more important at 2 /xm. The spectrum is indeed charac- 
terised by strong CO absorption longward of 2.3 /xm, indi- 
cating that the effect described above is quite significant at 
that wavelength. 

Another important consideration is whether the neb- 
ulae are ionisation bounded or density bounded. The for- 
mer assumption simplifies the analysis of emission lines, but 
there is some evidence that a significant number of the UV 
photons could escape the star-forming cloud without being 
absorbed, as has been suggested for NGC 4214 by Leitherer 
et al. (1996). In such cases the size of the Hll region is limited 
not by the number of ionising photons but by the amount 
of gas available which can be ionised, so although there may 
be O stars producing a large Lyman continuum (Lyc) flux, 
the observed Balmer line fluxes are all less by a constant fac- 
tor (leaving line ratios unaltered). On the other hand, since 
the ionisation potential for O+ of 35.1 eV is higher than the 
13.6 eV for H°, the [Om] lines originate from much closer 
to the star cluster, and so are unlikely to be truncated in 
the same way as the hydrogen lines. Thus for the same Lyc 


flux, a density bounded nebula would have a higher ratio of 
[Om]/Ha than an ionisation bounded nebula. Observation 
of a low ratio strengthens the evidence for a lack of O stars. 
The UV continuum leads us to a similar conclusion since, 
as discussed in Section 3, the observed continuum has the 
same shape as the intrinsic continuum which indicates an 
origin from A or B type stars. 

The third possibility we consider is that dust mixed 
within the nebulae is competing for UV photons, another 
effect which is not detectable from Balmer line ratios. Al- 
though this may be more significant than extinction beyond 
the limits of the nebula, it is often ignored and may be why 
models sometimes over-estimate the Balmer line fluxes, par- 
ticularly in young starbursts. Wood & Churchwell (1989) 
compared the radio contiuuum and far-infared fluxes in a 
number of ultra-compact Hll regions, and concluded that 
50-90 per cent of the Lyman continuum photons were ab- 
sorbed by dust internal to the nebulae. It may be that the 
discrepancy between the estimated total Ha flux and the 
measured surface brightness can be resolved if 80-90 per 
cent of the Lyc photons are absorbed by dust internal to 
Hll regions in NGC 7771. Dust absorption of Lyc photons 
is important even at optical extinctions low enough that 
relatively few [Om] and Ha photons are absorbed. Since 
the [Om] line originates in only the inner most part of the 
Hll region where the temperature is probably high enough 
to evaporate dust grains, its flux will be only marginally af- 
fected by internal dust. On the other hand the Ha line is 
produced throughout the nebula and will be attenuated, by 
absorption of Lyc photons. Thus for the same emitted Lyc 
flux, a dusty nebula will have a higher [OlIl]/Ha than a neb- 
ula with no internal dust. Again, observation of a low ratio 
strengthens the evidence for a lack of O stars. 

Next we construct a synthetic spectrum using a stellar 
library to match the starburst in Region 1, after correction 
for the evolved stellar population in the galaxy ’s bulge. To 
do this we extracted a third spectrum from the region 5.25- 
9.75 arcsec northeast of the centre which we showed to have 
an old red population (see Section 2.2). After correcting for 
the different aperture sizes, this spectrum was subtracted 
from that of Region 1 assuming that they both suffer the 
same extinction of Av = T9 mag. The result was compared 
to the synthetic spectrum after it had also been reddened 
by Av = 1.9 mag, as shown in Figure 9. We find that the 
population of the starburst is consistent with domination by 
B stars. Any attempt to include O stars in the simulation 
would then also require the inclusion of K or M stars to 
retain the slope of the continuum, but the spectrum shows 
little evidence for absorption features associated with such 
stars. The scenario of a population dominated by B type 
rather than O type stars, is consistent with the low [Om] 
flux and [Om]/Ha ratio, as well as the slope of the UV 
continuum. 

A simple explanation for the dearth of O-stars is that 
starburst is older than a typical O-star, and that starforma- 
tion has now ceased. Models such as the one described can 
be used to estimate the time evolution of EWh q , including 
contributions to the continuum from the nebular emission 
as well as the starburst population. Alonso- Her rero et al. 
(1996) include a contribution from an evolved underlying 
stellar population, but as we have discussed this is not sig- 
nificant at optical wavelengths for NGC 7771. We find that 
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Figure 9. Solid line: spectrum of the starburst (Region 1) with the evolved bulge population subtracted. The short wavelength end 
(< 4000 A) in particular is dominated by the noise from the latter. Dotted line: synthetic spectrum of 3 types of B star from the library 
of Silva & Cornell (1992) with Ay =1.9 magnitudes of extinction. Contributions to the flux at 5450 A are 60 per cent B5III, 20 per cent 
Bll t and 20 per cent B8l. 


in order to predict EWh q ^ 100 A the starburst must have 
ceased 6-8 million years previously. 

One problem with this idea is that the optical morphol- 
ogy suggests that it is a merger in its early stages (Sec- 
tion 5.4), typical of systems in which star-formation is just 
beginning rather than ending. The barred gravitational po- 
tential acts so as to drive gas inwards by removing angu- 
lar momentum. The resulting increase in the cloud collision 
rate, which is proportional to the square of the cloud num- 
ber density, has two efFects: giant molecular clouds are more 
likely to be formed by coalescence, and compression induced 
cloud collapse is enhanced (see for example Telesco et al. 
1991, and a review of the emipirical law SFR oc with 
1 < n < 2 in Kennicutt 1989). The result of this widespread 
density increase is that star-formation occurs over very ex- 
tended scales. 

What process might stop this star-formation? Clearly, 
if the gas supply ceases, so will the star-formation, but 
NGC 7771 has gas in abundance. The mass of molecular 
hydrogen is estimated from the CO luminosity (measured 
with a 55 arcsec beam) to be log (Mh 7 /Mq) = 10.3 (Sanders 
et al. 1991). An alternative might be a superwind created 
by the piling up of efficiently thermalised supernovae which 
could ionise the ISM, and perhaps the ensuing mass outflow 
could halt the inflow. Superwinds are discussed in detail by 
Heckman et ad. (1990) who find outflows along a galaxy’s mi- 
nor axis, and therefore this is unlikely to hinder the inflow 
of fresh gas in the disk plane. We nevertheless examine the 
evidence for a superwind in NGC 7771. Note: unfortunately 
the orientation of the slit is not ideal for observing a wind 
along the galaxy’s minor axis. One piece of evidence would 
be optical emission lines with ratios indicating shock-heated 
gas, in particular strong [Sll] and [Ol] lines. Table 2 gives 
[Sll](A6716 + A6731)/Ha = 0.18 and [Ol]A6300/Ha = 0.02, 
rather less than typical shock-excited values of 0.03-0.50 and 
0.25-1.20 respectively. Due to the inclination of NGC 7771 
we might also expect to see a kinematic signature of double 
peaked line profiles with a line splitting of 200-600 km s" 1 , 
this is not observed. From the correlations in Heckman et 
al. (1990) between line strength and width we would predict 
line widths around 50 km s -1 , consistent with some of our 


observations, but implying very little mass outflow. Finally 
we recall that the Ha flux is associated with a strong blue 
continuum, which would not be the case if the line emission 
has an origin in shock heated gas rather than stars. Thus, at 
least in the central few kiloparsecs, a superwind is not dy- 
namically or radiatively important. We also consider X-ray 
emission, which would be produced in a wind-blown cavity 
that has been heated to around 10 7 K. Although the reso- 
lution of the X-ray image is limited, it shows evidence for 
extended emission along the disk (as does the radio contin- 
uum), indicating star- formation within the bar, rather than 
perpendicular to the disk on the scale of tens of kiloparsecs 
as seen in some outflows. Thus we can say that a wind, if 
present at all, does not dominate the kinematics and lumi- 
nosity of NGC 7771. 

We therefore consider that unusual conditions would be 
required to halt the star- formation in its early stages, and 
conclude that it is still in progress. The corollary of this is 
that the starburst cannot be forming very massive stars. 


5.3 Truncated Mass Function 

From the evidence we have discussed in the previous section, 
we conclude that there are far fewer O stars than would 
normally be found in a starburst. This conclusion is based 
on several lines of argument: (1) the low equivalent width of 
Ha which cannot be explained in terms of geometry, density 
bounded nebulae, or dust absorption within the nebulae; (2) 
the weak [Om] lines, implying that the UV radiation field is 
not hot enough to ionise much 0 + ; (3) the UV continuum 
being very weak and having a slope typical of A or B stars; 
(4) the optical continuum of the star-forming region being 
dominated by B stars. Further we find no evidence that the 
star-formation should have ceased, and are led to conclude 
that no O stars are being made. 

We have re-calculated the starburst parameters for a 
Salpeter IMF but using an upper mass cutoff of 30 Mq for 
the model described previously and find that with continu- 
ous star-formation we can reproduce EWh q « 100 A. 

If we again assume constant star- formation on 
timescales of 10 8 yrs then r/ SN ~ 0.8 yr“ J , in much bet- 
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ter agreement with the radio and X-ray estimates than be- 
fore. The lower mass limit of 1 M© used in the models gives 
SF R ~ 40 Mq yr” 1 ; the effect of having a lower mass limit 
of 0.1 Mq would be to increase the mass of stars formed 
without significantly affecting the luminosity. We can inves- 
tigate the importance of this by considering the total mass of 
stars formed, assuming that stars greater than 8 Mq return 
all their mass to the interstellar medium (ISM) via winds 
and supernovae, while stars less massive than this effectively 
have all their mass removed from the ISM permanently. A 
lower mass limit of Afi = 1 M© would mean that 75 per cent 
of the burst mass is retained in stars, leading to a depletion 
by about 15 per cent of the current gas mass. Similarly we 
find that for M\ =0.1 Af©, as much as 90 per cent of the 
burst mass is removed permanently as it is tied up in low 
mass stars, but still leading to a depletion of only 20 per 
cent of the gas mass over the same period. The large re- 
serves of molecular hydrogen available suggest that either of 
these scenarios is equally feasible. 

The total Ha flux should be less than before since the 
most massive stars that dominated the ionising flux are 
no longer included in the calculation. We estimate it as 

4.3 x 10 -15 W m” 2 , which would extend over an area of 
60 arcsec 2 for our measurement of the Ha surface bright- 
ness. If we assume that the emission is extended by similar 
amounts both along and perpendicular to the slit, then this 
measurement is completely consistent with the observed ex- 
tent of 7.5 arcsec. 

5.4 Interaction 

NGC 7771 is part of an interacting system, being closely 
associated with NGC 7770 which is a smaller face-on galaxy. 
Their nuclei are separated by only 65 arcsec, a projected 
distance of 27 kpc, and their redshifts differ by ~350 km s -1 . 
Veilleux et al. (1995) have presented long-slit data on this 
galaxy (denoted by them NGC 7771 S) which show it has 
a similar Ha surface brightness, 27 x 10” 17 W m -2 in a 
2 arcsec x 2 arcsec aperture and corrected for Av ~ 1.8 mag 
extinction. The extent of this starburst is unknown, but the 
radio continuum shows the total flux is only about 10 per 
cent of that in NGC 7771. 

Mihos & Hernquist (1994) have modelled ‘minor merger 
systems’ which involve a disk galaxy accreting a low-mass 
companion. The tidal perturbations induced by the satel- 
lite can cause the disk galaxy to develop a two-armed spiral 
pattern which drives a substantial fraction of the disk gas 
towards the centred regions, triggering star formation. The 
presence of a bulge in the disk galaxy will help to stabilise 
the disk against strong inflow so that the final nuclear star- 
burst is less compact. 

Qualitatively, this model describes the situation in 
NGC 7771 reasonably well. It is interacting with a smaller 
galaxy, and it does appear to have a two-armed spiral; there 
is star-formation occuring weakly throughout the galaxy and 
in a much greater concentration towards the centre; the star- 
burst is certainly not compact, covering as it does a region 
more than 3 kpc across. 

We now attempt to explain the existence of the clus- 
ter of non-ionising stars that are obscuring part of the nu- 
clear starburst. This could be due to resonances associ- 
ated with a barred gravitational potential. If we consider 


the epicyclic approximation for stellar orbits, where a ret- 
rograde elliptical motion is superimposed on a prograde cir- 
cular orbit, the stars at some radius r will have an angular 
velocity Q(r) perturbed by an epicyclic frequency «(r) where 
«(r) 2 = -|- 40 2 (see Binney &: Tremaine 1987, Combes 

et al. 1995); we denote the pattern speed of the slowly rotat- 
ing bar to be f !*>. If inner Lindblad resonances were to exist 
they would do so where Qb = Q(r)-x(r)/ 2, although as this 
appears not to be so then at all radii Qb > Q(r) - k(t)/ 2 . 
The stars will also be oscillating perpendicular to the plane 
according to z = — v\z where the frequency u z is related the 
gravitational potential U(r,z) by v\ = BU ^ 0 ^ (Combes et 
al. 1995). There will be resonances, analogous to the Lind- 
blad resonances, when stars cross the disk plane in approx- 
imately the same position relative to the bar each time, ie 
Qb = H(r) — Vz{ r )/2. This can be at a different radius to 
the Lindblad resonances, and can be satisified even if the 
criterion for Lindblad resonances is not. It has been used to 
explain the peanut-shape of some edge-on galaxies, and in 
principle it could be that the obscuring cluster is a deforma- 
tion due to such a resonance. The scenario would account 
for the higher velocity dispersion for the emission lines in 
this region. The typical timescale for the stars to travel a 
distance on the order of 1 kpc above the disk plane (similar 
to their distance from the centre in the plane), with a typical 
velocity of 100 km s“ 1 (from Section 2.3) is 10 7 years, suffi- 
cient for stars more massive than ~ 10 M© to have evolved 
away from the Main Sequence. Thus we would expect, as is 
observed, that the obscuring stars do not include any with 
strong ionising fluxes. 


6 CONCLUSION 

We present new optical and UV spectra, and a soft X-ray 
image and spectrum, which are analysed in conjunction with 
data taken from the literature, covering the entire spectral 
energy distribution of the galaxy. 

The optical spectrum has shown that there is a large 
starburst in the nucleus of NGC 7771 that is partially hid- 
den by a cluster of non-ionising stars. Although the infrared 
(8-1000 ^m) and soft X-ray emission are confused by the 
proximity of NGC 7770, the radio flux shows that this lat- 
ter galaxy makes only a small contribution. Thus we can 
adopt a supernova rate of 0.8-1 yr -3 as indicated by these 
three fluxes. The estimated total Ha flux is consistent with 
the measured Ha surface brightness and its extent along 
the slit, assuming that the emission is extended a similar 
amount perpendicular to the slit. 

The low equivalent width of Ha, on the order of 100 A, 
is apparently at odds with the high current star-formation 
rate. By also considering other evidence — namely the weak 
[Om] lines, the shape and weakness of the UV continuum, 
and the shape of the optical continuum in the starburst — 
we have argued that this cannot be due to effects of geome- 
try, density bounded nebulae or dust absorption within the 
nebulae. Modelling the starburst with a standard IMF indi- 
cates that the star formation episode must have concluded 
several million years previously, a result that is inconsistent 
with the large mass of available gas and the dynamical state 
of the galaxy. We conclude that the initial mass function is 
truncated, with M u = 30Af©. 
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NGC 7771 is in the early stages of a merger, the massive 
burst of star formation currently in progress apparently hav- 
ing been caused by tidal interaction with the nearby dwarf 
galaxy NGC 7770. This interaction has given rise to a barred 
potential which is driving gas inwards and fuelling the star- 
burst. The bulge in NGC 7771 acts as a regulator, opposing 
the inflow caused by the interaction and hence the star- 
formation is extended over several kiloparsecs. This may also 
cause the anomalously low fraction of massive stars com- 
pared to the standard picture of a starburst. 

Near the centre there is a resonance, analogous to Lind- 
blad resonances but perpendicular to the disk plane. This is 
causing stars in the region to oscillate out of the plane more 
vigorously and gives rise to what appears to be a cluster of 
stars obscuring part of the starburst. The timescale for the 
stars to move out of the plane acounts for the lack of ionising 
stars in this cluster. 
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